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Uncertainty  in  Null  Polarimeter  Measurements 

K.  B.  Rochford  and  C.  M.  Wang 

National  Institute  of  Standards  and  Technology 
325  Broadway,  Boulder,  CO  80303-3328 


This  Internal  Report  details  the  error  analysis  of  a nulhng  polarimeter  used  for  retardance 
measurements.  We  determined  that  the  uncertainty  arising  from  random  effects  is 
between  0.07°  and  0.10°  for  measurements  of  several  retarders  with  nominally  90° 
retardance.  This  instrument  and  two  other  methods  were  used  to  determine  the  retardance 
of  a stable  rhomb  device  proposed  as  a standard  retarder.  The  measurement  results  were 
used  to  support  certification  of  our  retarder  as  a NIST  Standard  Reference  Material. 

Keywords:  Fresnel  rhomb;  measurement;  null  polarimetry;  polarimetry;  polarization; 
retardance;  standards;  Standard  Reference  Materials;  uncertainty  analysis. 


1.  Introduction 

A system  commonly  known  as  the  manual  null  instrument  [11  or  nulling  elhpsometer  [2]  was  a secondary 
polarimetric  method  used  for  determining  retardance  in  our  efforts  to  develop  a stable  retarder  [3,4].  For 
this  measurement,  a compensating  quarterwave  retarder  (C)  and  the  retarder  under  test  (R)  are  placed 
between  an  input  polarizer  (P)  and  analyzing  polarizer  (A)  (Figure  1).  We  illuminate  the  system  with 
circular  polarization  so  that  the 
optical  power  incident  on  the 
compensator  is  constant  as  the 
azimuth  of  P is  rotated.  This  circular 
input  polarization  is  obtained  using  a 
linearly  polarized  (LP)  optical  beam 
and  a Soleil-Babinet  compensator 
(SBC)  adjusted  to  quarterwave 
retardance. 

To  measure  retardance,  the  system  is  aligned  so  that  polarizers  P and  A are  crossed  (transmission  axes  are 
orthogonal)  and  the  retardance  axes  of  compensator  C are  coincident  with  the  polarizer  axes.  The 
retardance  axes  of  R are  aligned  at  45°  to  the  polarizer  axes.  While  the  optical  power  exiting  the  system  is 
monitored,  polarizer  P is  rotated  to  the  angle  that  yields  the  minimum  detected  power.  The  magnitude  of 
polarizer  rotation  a required  to  reach  this  null  is  equal  to  half  the  measured  retardance  6^  (or  6^  = 2a). 

This  system  benefits  from  the  high  accuracy  with  which  transmission  nulls  can  be  determined.  If  the 
system  were  used  without  the  compensator  C,  the  output  of  retarder  R would  be  elliptical,  in  general,  and 
the  null  could  be  difficult  to  locate.  For  example,  if  R is  a quarterwave  retarder,  the  system’s  transmittance 
is  constant  for  all  positions  of  P if  no  compensator  is  used,  and  a null  cannot  be  found.  Because  an 
appropriately  ahgned  quarterwave  retarder  can  convert  any  elhptical  state  into  a linear  state,  component  C 
causes  the  hght  incident  on  the  polarizing  analyzer  to  be  linearly  polarized,  and  this  ensures  a high- 
extinction  null  that  can  be  accurately  determined. 
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Figure  1 . Schematic  of  null  polarimeter  apparatus. 
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Generally,  several  parameters  are  important  for  a complete  ellipsometric  measurement.  For  example,  the 
sign  of  the  compensator’s  retardance  (location  of  the  fast  axis)  and  the  direction  of  polarizer’s  rotation  are 
required  to  determine  the  sign  of  the  retarder.  In  addition,  the  analyzer  can  be  rotated  to  a null  and  provide 
a measure  of  diattenuation.  For  our  purposes,  however,  these  details  were  unimportant  and  ignored.  First, 
the  sign  of  the  retardance  is  determined  by  the  design  of  our  rhomb  retarder,  and  routine  experimental 
determination  is  not  needed.  Second,  the  device  is  designed  for  negUgible  diattenuation,  and  this  has  been 
verified  independently,  so  additional  measurement  is  not  necessary. 

In  this  paper,  we  discuss  the  retardance  errors  associated  with  the  null  measurements  made  on  five 
prototype  stable  retarders  labeled  SRI  through  SR5. 

2.  Experimental  Method 

Each  retardance  measurement  required  several  alignment,  readout,  and  adjustment  steps.  The  procedure 
used  to  obtain  a retardance  estimate  6^,  is  discussed  below.  Also,  several  general  concepts  and  practices 
are  reviewed. 

2.1  General  Practices 

2.1.1  As  optical  elements  are  inserted  into  the  optical  beam  path,  the  elements  are  adjusted  for  normal 
incidence  by  monitoring  the  Fresnel  reflection  from  the  surface.  For  anti-reflection  coated 
components,  this  was  done  using  a colhnear  HeNe  laser  directed  along  the  IR  laser  path. 

2.1.2  “Nulling,”  or  finding  the  polarizer/analyzer  angle  which  yields  maximum  extinction,  is  performed 
several  times  over  the  course  of  a measurement.  Two  methods,  creatively  named  “A”  and  “B,” 
are  described  below; 

Method  A:  The  device  (polarizer  P or  analyzer  A)  is  manually  rotated  to  the  point  of 
maximum  extinction.  The  angle  is  read  directly  from  the  vernier  on  the  rotation  stage. 

Method  B:  The  device  is  rotated  to  a position  that  yields  a predetermined  transmittance  (for 
example,  T = lO"*),  and  the  angle  read  directly  from  the  stage  and  recorded.  The  device  is  then 
rotated  through  the  null  to  a second  point  with  the  same  transmittance,  and  this  second  angle 
recorded.  This  procedure  may  be  repeated  several  times.  The  arithmetic  mean  is  calculated, 
and  the  stage  is  rotated  to  the  mean  angle. 

2.1.3  The  procedure  requires  that  the  retardance  axes  of  several  components  be  found  polarimetrically. 
When  a retarder  is  placed  between  a crossed  polarizer  and  analyzer,  minimum  transmission  occurs 
when  the  retardance  axes  coincide  with  the  polarizer/analyzer  axes.  Thus,  rotating  the  retarder 
until  the  system  output  is  nulled  locates  the  retardance  axes. 

2. 1 .4  The  input  to  the  first  polarizer  is  collimated,  circularly  polarized  hght  from  a laser  diode. 
Collimated  light  was  polarized  and  directed  toward  an  analyzing  polarizer  (A),  lens  and  linear 
photodetector  aligned  with  the  propagation  axis  (the  Figure  1 setup  with  components  P,  C,  and  R 
removed  from  beam  path).  SBC  is  placed  after  the  first  polarizer,  and  its  retardance  axis  located 
polarimetrically.  The  SBC’s  retardance  axes  are  rotated  45°,  and  the  retardance  of  the  SBC  is 
adjusted  to  quarterwave.  Proper  SBC  adjustment  is  checked  by  measuring  the  modulation  of 
detected  power  as  the  analyzer  is  rotated;  there  should  be  no  modulation  for  perfect  circular 
polarization.  Once  aligned,  the  optics  for  providing  circular  polarization  were  not  routinely 
adjusted  between  retardance  measurements. 
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2.2  Experimental  Procedure 

The  procedure  for  measuring  retardance  consists  of  the  seven  steps  described  below; 

2.2.1  A polarizer  (P)  is  inserted  into  the  system  after  the  SBC.  The  transmission  axis  is  aligned  roughly 
perpendicular  to  the  optical  table,  and  the  analyzer  (A)  is  rotated  to  a transmission  null  (maximum 
extinction).  The  polarizer  angle  is  read  to  obtain  the  initial  angle  ajo.  Nulling  and  reading  are 
accomplished  using  either  method  A or  B,  which  are  used  throughout  the  procedure  and  are 
described  above  (Section  2. 1 .2). 

2.2.2  The  retarder  R is  inserted  into  the  system  and  its  retardance  axes  found  polarimetrically.  Then  the 
retarder  is  temporarily  removed  from  the  system.  (The  mounting  hardware  is  sufficiently  stable  to 
allow  extraction  and  reinsertion  without  significantly  altering  retarder  alignment.) 

2.2.3  The  polarizer- analyzer  pair  are  rotated  45°  so  that  the  polarizer- analyzer  transmission  axes  bisect 
the  retardance  axes  of  R.  The  initial  angle  for  retardance  measurement  becomes  a,  = Kjo  + 45°. 

2.2.4  The  compensating  quarterwave  plate  (C)  is  inserted  near  polarizer  P.  The  compensator’s 
retardance  axes  are  aligned  to  coincide  with  the  polarizer- analyzer  axes. 

2.2.5  The  retarder  is  reinserted  in  between  compensator  C and  analyzer  A.  Polarizer  P is  rotated  to  a 
transmission  null  which  occurs  at  polarizer  angle  ttj.  The  measured  retardance  6^,,  = 2 1 I • 

2.2.6  Components  C and  A are  rotated  90°  so  that  the  measurement  can  be  repeated.  (This  second 
measurement,  when  averaged  with  the  results  of  the  prior  measurement,  cancels  the  effect  of 
optical  activity  in  compensator  C and  has  no  detrimental  effect  on  measurement  accuracy.) 

2.2.7  Polarizer  P is  rotated  to  a transmission  null,  and  another  a2  is  found.  A second  measured 
retardance  6^^  = 2|  a,  - a,  | is  calculated.  The  measured  retardance  of  our  specimen  R is  6^ 

=(6«+6„2y2. 

2.3  Experimental  Limitations 

This  procedure  uses  real  components  that  are  imperfect  or  have  limitations.  Several  of  the  component 

issues  of  interest  follow: 

2.3.1  The  optical  source  is  a packaged  laser  diode.  The  device  is  driven  with  constant  current,  but  still 
undergoes  changes  in  wavelength  due  to  temperature,  age,  etc.  We  estimate  the  worst-case 
wavelength  variation  to  be  1.31  ± 0.02  pm. 

2.3.2  The  input  polarization  is  not  perfectly  circular.  When  light  was  incident  upon  a rotating  polarizer, 
a worst-case  modulation  of  20%  could  be  measured.  This  corresponds  to  a SBC  retardance  of  90 
± 6.4  °.  (SBC  retardance  653  = cos“’[(In,^  - L.n)/(Imax  + Imin)  = cos’ '[(1.0  - 0.8)/(1.0  -1-  0.-8)] ). 

2.3.3  The  polarizers  used  in  the  systems  are  high  quality  calcite  polarizers.  For  our  measurements,  the 
light  leakage  at  extinction  overwhelmed  any  detection  noise,  and  the  polarizers  limit  the  null 
determination.  The  maximum  extinction  was  found  by  measuring  the  minimum  and  maximum 
transmission  through  crossed  polarizers  with  phase-sensitive  detection.  Four  measurements  of  the 
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extinction  (=  Imin/Imax)  for  the  polarizers  P and  A yielded  a mean  of  6. 1 x lO'^  with  a standard 
deviation  of  1.9  x 10'^. 


2.3.4  The  retardance  of  compensating  waveplate  C was  independently  measured  to  be  87.9°.  Including 
the  effect  of  ± 20  nm  wavelength  variation,  the  range  of  the  compensator’s  retardance  is  87.9 


± 1.4°. 


2.4.5  The  rotation  stages  used  in  the  experiment  were  ruled  with  marks  at  each  0.1°.  The  rotation  was 
read  directly  from  etched  hnes  on  the  rotating  and  fixed  mounts  (and  not  read  through  any  geared 
mechanism  or  screw  drive),  so  the  measurements  were  free  of  backlash. 

3.  Mathematical  Modeling 

For  optical  systems  that  have  no  depolarizing  components,  the  evolution  of  polarization  can  be 
satisfactorily  modeled  using  Jones  calculus  [5].  The  measurements  of  our  rhomb  retarders  in  the  null 
polarimeter  meet  this  requirement,  and  Jones  calculus  was  used  to  model  the  ideal  system  and  determine 
the  uncertainty  due  to  individual  component  imperfections  and  ahgnment  uncertainty. 

In  Jones  calculus,  the  change  in  the  electric  field,  defined  as  the  column  vector  [E^^,  Ey]^,  is  found  by 
multiplying  the  input  field  vector  by  a 2x2  matrix  corresponding  to  an  optical  element.  Thus  the  output 
vector  Eq  resulting  from  a field  Ej  propagating  through  optical  elements  X and  Y (which  have  Jones 
matrices  M,;  and  My)  is  Eo  = My  x M^  x E;. 

3.1  Ideal  Null  System 

For  the  ideal  null  system,  the  following  Jones  matrices  and  vectors  adequately  describe  the  components  of 
interest: 

3.1.1  Light  entering  the  measurement  system  is  circularly  polarized  with  input  field 


(1) 


3.1.2  Polarizer  P at  angle  a from  origin 


cos^(a)  sin(a)cos(a) 

P = 

sin(a)cos(a)  sin^(a) 


(2) 


3.1.3  Compensator  C with  90  ° retardance 


C = 


(3) 


0 e 


3.1.4  Retarder  R rotated  45°  with  unknown  retardance 
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RET- 


(4) 


cos(6q/2)  /sin(6Q/2) 
/sin(6Q/2)  cos(6q/2) 


3.1.5  Analyzer  A fixed  at  90° 


0 

1 


The  output  field  is  found  by  matrix  multiplication 

E^=AxRETycCxP^E. 


(5) 


(6) 


and  the  output  intensity  is  calculated  from  the  vector  Eq  using  I = E^E^’  + EyEy*,  where  the  asterisk  denotes 
complex  conjugation.  For  the  ideal  system  described  by  the  matrices  above, 

1 -(sin(2a)sin(6Q)  +cos(2cc)cos(6q)  ) 

i 


The  null  occurs  where  dl/da  = 0,  or  when  tan(2a)  = tan(6o);  thus,  = 2a  for  the  ideal  system. 

3.2  Measurement  error  due  to  individual  component  error 

Several  possible  error  sources,  including  component  imperfections  or  alignment  uncertainty,  were 
identified  after  careful  consideration  of  the  null  system.  Each  of  these  errors  is  analyzed  by  forming  the 
Jones  matrix  corresponding  to  the  component-  alignment  error  and  determining  the  effect  on  measured 
retardance.  Six  errors  were  identified  and  analyzed  below: 

3.2. 1 Error  due  to  offset  in  compensator  retardance.  If  the  compensator  retardance  is  nil  + €,  the  Jones 
matrix  becomes 


C = 


gi(7r/2^e)/2 

0 


0 


g -;(7t/2+e)/2 


(8) 


and  the  detected  intensity  as  a function  of  polarizer  position  becomes  I = (1  - sin(2a)sin(6o)cos(e) 
- cos(2a)cos(6o))/4.  Again  the  null  occurs  when  the  derivative  dl/da  = 0,  or  when  tan(2a) 

= tan(6o)cos(€).  Since  the  measured  retardance  6^,  = 2a,  tanCb^j)  = tan(6Q)cos(e),  and  the 
retardance  measurement  error  A6  = bg  - 6^,  = bg  " tan'‘[tan(bo)cos(6)]. 
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3.2.2  Error  due  to  compensator  misalignment.  If  the  compensator  retardance  axes  and 

polarizer- analyzer  axes  are  misaligned  by  e,  the  Jones  matrix  for  the  compensator  becomes 


e''^'^cos\e) +e  /sin(2€)sin(Tt/4) 

isin(2e)sin(7t/4)  e +e'’^^'^sin^(€) 


(9) 


Again,  we  calculate  the  detected  intensity  I and  find  where  the  derivative  dl/da  = 0.  This  yields 

cos(6Q)sin(4€)  +2sin(6Q)cos(26) 


tan(2a)  = 


cos(6 J cos(4€)  - 2 sin(6 J sin(2e)  + cos(6 J 


(10) 


As  before,  the  error  in  measured  retardance  A 6 can  be  found  by  subtracting  2a  ( = 6^)  from  Sq. 

We  have  treated  both  retardance  and  alignment  compensator  errors  as  separate  here  and  in  Section 
3.2.1.  When  both  errors  are  taken  together, 

cos(6)sin(4e.,)[  1 +sin(e,)]  +2sin(6)cos(e,  )cos(2€,) 

tan(2a)= , ni) 

cos(6)[cos(4e2)[sin(e,)  + 1]  -sin(ej)  +1]  -2sin(6)sin(2e2)cos(ej ) 


where  €[  = is  the  compensator  retardance  error  and  €2  is  the  alignment  error. 

For  the  errors  e, , typical  for  our  measurements,  the  retardance  error  found  using  eq  (1 1)  is  equal 
to  the  sum  of  the  separate  errors  (eq  (10)  and  Section  3.2.1).  We  use  the  separate  forms  as  this 
allows  us  to  use  different  types  of  errors  e in  the  calculation  of  uncertainty. 


3.2.3  Error  due  to  misaligned  retarder.  If  the  retarder  is  misahgned  by  an  angle  e,  the  appropriate  Jones 
matrix  is 


RET  = 


'^‘^"sin‘(e  + Tt/4)  +e'^'^"cos^(e  +Tt/4) 
/sin(6Q/2)cos(2e) 


tsin(6Q/2)cos(2€) 

e '^‘^^sin^(e  + 7i/4)  +e  ‘^’^^cos\e  + Tt/4) 


(12) 


Calculating  the  derivative  dl/da  gives 

2sin(6jcos(2€) 

tan(2a)  = . 

cos(4e)(cos(6Q)  - 1)  +cos(6q)  + 1 


(13) 


Measured  retardance  error  A6  can  be  found  by  subtracting  2a  ( = 6n,)  from  6q  This  error  causes  a 
decrease  in  measured  retardance  for  all  cases  of  interest. 
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3.2.4 


Error  due  to  analyzer  misalignment.  If  the  analyzer  is  misaligned  by  an  angle  e,  the  Jones  matrix 
becomes 


A = 


sin^(e)  -cos(6)sin(e) 


(14) 


-cos(e)sin(€)  cos^(€) 


Solving  as  before,  we  find  that  the  dependence  on  e cancels  and  that  tan(2a)  = tan(6),  the 
relationship  for  an  ideal  system.  Thus,  errors  in  analyzer  position  will  not  affect  retardance 
measurements. 

3.2.5  Error  due  to  input  polarization  error.  If  the  optical  input  to  the  system  is  not  circularly  polarized, 
errors  can  result  when  the  null  is  found  using  method  B.  This  occurs  because  the  transmission 
through  polarizer  P does  not  vary  symmetrically  with  polarizer  angle  if  the  input  polarization  is  not 
circular. 

This  can  be  modeled  by  using  an  input  field  matrix  that  is  a function  of  SBC  retardance  653 


(15) 


To  find  the  functional  form  of  intensity  as  the  polarizer  rotates,  we  model  the  simpler  system  I = A 
X P X Ej,  which  yields  I = cos(a)sin^(a)cos(6sB)  + l/2sin^(a).  For  Inu,,  (which  corresponds  to  the 
intensity  at  which  angles  on  either  side  of  the  null  are  measured)  and  retardance  653  , we  can  find 
the  roots  a'  and  a".  Comparing  the  average  = (a'+  a")/2  with  the  ideal  value  a = 6(/2  gives 
the  angular  error  caused  by  input  polarization  errors. 

3.2.6  Error  due  to  optical  activity  in  compensator  waveplate.  The  quarterwave  compensator  had  a small 
amount  («  1 °)  of  optical  activity  (or  rotation).  This  was  modeled  using  a Jones  matrix  for  C that 
included  both  retardance  and  rotation.  Though  the  analysis  shows  that  significant  error  can  occur 
when  determining  the  retardance  (Section  2.2.1  through  2.2.5),  the  sign  of  the  error  changes  when 
the  retardance  is  remeasured  after  rotating  the  compensator  and  analyzer  by  90°  (Section  2.2.6). 
Thus,  averaging  the  two  retardance  measurements  cancels  this  error,  and  the  null  measurement  is 
insensitive  to  optical  activity  in  the  compensator. 

4.  Analysis  of  Measurement  Uncertainties 

We  can  calculate  retardance  uncertainty  for  each  step  of  a retardance  measurement.  As  slightly  different 
component  values  were  used  for  individual  measurements,  retardance  uncertainty  for  individual  rhombs  is 
hsted  as  required. 

4.1  General  uncertainties 

The  measurement  of  retardance  requires  the  measurement  of  rotation  angles  as  well  as  accurate  rotation  of 
a components  to  specified  angles.  We  denote  the  calculated  angle  uncertainty  by  e^,  where  the  subscript  x 
labels  the  source  of  error. 
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4.1.1 


Stage  resolution.  The  rotation  stages  used  for  these  measurements  were  ruled  every  0.1  ° and  were 
read  from  the  overlap  of  graduations  on  two  ruled  pieces  (one  fixed  vernier,  and  one  rotating 
scale).  We  assume  that  during  readout  we  can  resolve  reading  within  ± 0.05°,  or  halfway  between 
the  graduations.  The  uncertainty  is  estimated  by  assuming  that  the  readout  has  a uniform  error 
over  of  a width  of  0.05°.  With  this  model,  the  standard  angle  uncertainty  (one  standard  deviation) 
is  6R  = 0.05°/yi2  = 0.0144°. 

4. 1 .2  Method  A nulling.  Transmission  nulling  using  method  A is  limited  by  polarizer  imperfections. 
Ideally,  the  transmittance  T = sin^O,  where  0 is  the  angle  from  extinction.  For  our  polarizers,  the 
best  extinction  obtained  was  T = 6.1x10'^  with  a standard  deviation  of  1.9x10'^  based  on  four 
measurements.  Thus  corresponds  to  an  angle  resolution  of  0.045°  (with  a standard  error  of 
0.0035).  We  calculate  the  maximum  angle  uncertainty  (3a)  as  the  mean  angle  resolution  plus  3 
standard  errors  (0.045°  + 3-0.0035°  = 0.0555°),  so  the  standard  uncertainty  for  nulling  = 
0.0185°. 

For  some  cases,  the  nulling  angle  must  be  read  out  for  calculation  or  angle  adjustment.  In  this 
case,  we  combine  the  nulling  and  readout  uncertainty  to  form  = 0.0235°. 

4.1.3  Method  B nulling.  Transmission  nulling  using  method  B is  more  complicated  because  of  input 
polarization  uncertainty  (which  is  due  to  uncertainty  in  the  SB  compensator  retardance).  If  the 
null  is  found  by  locating  the  angles  at  which  the  transmittance  is  Inui,  when  the  compensator 
retardance  is  653  (Section  3.2.5),  the  worst-case  uncertainty  Aesb  = (a'  -1-  a")/2  where  a'  and  a" 
are  the  roots  of  cos(a)sin^(a)cos(6sB)  -1-  ‘/2sin^(a)  - Iny,,  = 0.  The  angle  uncertainty  due  to  input 
polarization  is  655  = AesB/3.  Because  the  values  of  Inu,,  differed  among  the  various  measurements, 
its  values  and  corresponding  653  are  tabulated  in  Table  1 for  653  = 90  ± 6.4°. 


Table  1.  Angular  uncertainty  due  to  input  polarization  errors. 


Device 

p 

Maximum  Iny,, 

Roots  (rad) 

( ) 

esB  (°) 

SRI 

2.44  X 10'^ 

0.00698,  -0.00699 

2.86  X 10-^ 

9.55  X 10  ' 

SR2 

5.95  X 10  ' 

0.01090,  -0.01092 

5.73  X 10-^ 

1.91  X 10" 

SR3 

1.52  X 10-" 

0.01740,  -0.01747 

2.01  X 10-' 

6.68  X 10-" 

SR4 

1.05  X 10  ' 

0.04561,  -0.04608 

1.35  X 10  ' 

4.49  X 10  ' 

SR5 

4.40  X 10-' 

0.02957,  -0.02977 

5.73  X 10  ' 

1.91  X 10' 

In  addition  to  uncertainty  €53,  additional  uncertainty  results  from  the  readout  process;  the  two 
angular  positions  of  Inuj,  must  be  read,  and  the  stage  then  rotated  to  the  average  of  these  two  angles. 
These  uncertainties  can  be  combined  to  form  the  total  uncertainty  for  nulling  using  method  B,  €3 
= (€53^  -1-  363^2)'^.  Table  2 lists  the  resulting  uncertainty  for  each  device. 

Table  2. 

Angle  uncertainty  from  method  B nulling. 

SRI 

SR2  SR3 

SR4 

SR5 

0.0177 

0.0177  0.0177 

0.0182 

0.0178 
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4.2  Angle  Measurement  Uncertainty  Arising  from  Experimental  Limitations 

Each  retardance  measurement  has  several  steps  requiring  reading  and/or  adjustments  of  the  rotation  angles 
of  various  system  components.  In  this  section,  we  list  the  uncertainties  arising  at  each  step  in  order  of 
occurrence.  (For  reference,  the  analysis  is  numbered  to  correspond  with  the  procedural  steps  described  in 
Section  2.2.) 

4.2.1  Null  polarizer- analyzer  and  readout  angle  a,.  Both  nulling  methods  were  used  in  determining  a, 
for  individual  measurements.  Also,  one  determination  (SR4)  used  the  average  of  two 
measurements.  The  angle  uncertainty  associated  with  this  step  is  either  e^R  or  Eg,  depending  on 
the  nulling  method  used.  For  N measurements,  the  angle  uncertainty  Ej  = E^Rg/N’^.  The  methods 
and  results  are  tabulated  in  Table  3. 

Table  3.  Angle  uncertainty  in  polarizer- analyzer  readout. 


SRI  SR2  SR3  SR4  SR5 


Null  method 

A 

B 

A 

B 

B 

N 

1 

1 

1 

2 

1 

€1 

0.0235 

0.0177 

0.0235 

0.0129 

0.0178 

4.2.2  Location  of  the  retardance  axes  of  R.  Method  A was  used  for  each  determination.  Thus,  the  angle 
uncertainty  associated  with  axis  determination  is  Ej  = e^  = 0.0185°. 

4.2.3  45°  rotation  of  polarizer/analvzer.  This  step  requires  a rotation  of  45°  with  respect  to  the  retarder 
axes,  so  rotation  causes  an  effective  uncertainty  in  the  location  of  these  axes.  The  uncertainty  of 
retarder  axis  angle  E2  and  a polarizer/analyzer  readout  with  angular  uncertainty  Er  are  combined. 
The  combined  error  for  the  retarder  axes  after  the  polarizer  / analyzer  rotation  is  Ere^  = (e2^ 

+ Er^)’^  = 0.0235  for  all  measurements. 

In  addition,  the  angle  of  the  polarizer  with  respect  to  the  analyzer  accrues  uncertainty  during  the 
rotation  because  of  the  initial  nulling  limitation  and  readout  error  from  rotation.  The  combined 
angle  uncertainty  between  the  polarizer  and  analyzer  is  Ep/^  = + €r^)'^.  Table  4 lists  the 

resulting  uncertainties. 


Table  4.  Angle  uncertainty  due  to  polarizer/analyzer  misahgnment. 


SRI 

SR2 

SR3 



SR4 

SR5 

Ep/A  0.0275 

0.0228 

0.0275 

0.0194 

0.0229 

4.2.4  Location  of  compensating  waveplate’s  retardance  axes.  The  determination  of  the  waveplate’s  axes 
was  made  using  different  methods,  and  the  resulting  angle  uncertainty  E4  equals  e^  or  Eg, 
depending  on  the  nulhng  method  used.  In  one  case,  two  measurements  were  made  so  the 
uncertainty  is  reduced  by  2 ’^^.  The  resulting  angle  uncertainty  for  each  case  is  tabulated  in  Table 
5. 
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Table  5.  Angle  uncertainty  in  waveplate  axis  determination. 


SRI 

SR2 

SR3 

SR4 

SR5 

Null  method 

A 

A 

A 

B 

B 

N 

1 

1 

1 

1 

2 

^4 

0.0185 

0.0185 

0.0185 

0.0182 

0.0126 

4.2  5 Rotation  of  polarizer  to  null  to  find  and  cto  - Null  method  B was  used  for  each 

measurement  so  the  null  angle  uncertainty  = 6b  for  each  case.  Multiple  measurements  were  made 
in  several  cases,  so  the  effective  angle  uncertainty  65  = eg/N'^.  The  difference  - a.^  is  of 


interest,  and  the  combined  uncertainty  is 
uncertainties  for  finding  aj. 

= (65^  + 6,^)’'^. 

Table  6 lists  the  calculated  angle 

Table  6. 

Angle  uncertainty  for  polarizer  rotation  to  null. 

SRI 

SR2 

SR3 

SR4 

SR5 

N 

3 

3 

1 

3 

3 

£5 

0.0102 

0.0102 

0.0177 

0.0105 

0.0103 

0.0256 

0.0204 

0.0294 

0.0167 

0.0205 

4.2.6  Rotation  of  biasing  waveplate  to  repeat  measurement.  The  quarterwave  compensator  and  analyzer 
are  rotated  90°.  Compensator  rotation  introduces  readout  uncertainty  and  an  angle  uncertainty 
€4  due  to  the  initial  alignment  of  the  compensator  These  uncertainties  combine  to  form  an  angle 
uncertainty  for  the  compensator  position,  = (6^^  + 64^)'^.  The  uncertainty  in  analyzer  rotation 
does  not  contribute  to  retardance  uncertainty  (see  Section  3.2.4)  and  can  be  ignored.  The  values 
for  €5  are  tabulated  in  Table  7. 


Table  7.  Angle  uncertainty  in  compensator  position. 


SRI 

SR2 

SR3 

SR4 

SR5 

^6 

0.0235 

0.0235 

0.0235 

0.0233 

0.0191 

4.2.7  Remeasurement  of  nulling  angle  and  a^-  cL].  As  in  Section  4.2.5,  the  nulling  angles  are  all 
measured  by  method  B so  the  angle  uncertainty  is  63.  Multiple  measurements  were  made  in 
several  cases,  so  the  effective  angle  uncertainty  €7  = Cg/N'^.  The  difference  a,  - a,  is  of  interest, 
and  the  combined  uncertainty  is  6^^.  = (67^  + Table  8 hsts  the  calculated  angle 

uncertainties. 


Table  8. 

Angle  uncertainty  for  second  rotation  to  polarizer  null. 

SRI 

SR2 

SR3 

SR4 

SR5 

N 

2 

4 

1 

3 

3 

€7 

0.0125 

0.0088 

0.0177 

0.0105 

0.0103 

€ao' 

0.0266 

0.0198 

0.0294 

0.0167 

0.0205 
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4.3  Uncertainty  in  Measured  Retardance 

The  angle  uncertainties  and  component  value  errors  or  uncertainties  combine  to  yield  uncertainty  in  the 
measurement  of  retardance.  Some  of  the  uncertainties  depend  on  the  retardance  itself,  and  the  measured 
retardance  6^  is  used  for  this  estimate.  Because  the  uncertainties  are  generally  small,  the  uncertainty 
estimate  is  not  greatly  affected  by  using  the  measured  value  in  our  calculations  and  recursive  estimation  is 
not  needed.  The  standard  retardance  uncertainty  u ( = 6^,  - 2a,  where  2a  is  the  measured  retardance 
expected  from  the  error  equations,  and  Sq  = needed)  are  discussed  below. 

4.3.1  Retardance  uncertainty  due  to  retarder  / polarizer  misalignment.  The  retardance  uncertainty  Upj-j 
due  to  retarder  misalignment  is  found  using  eq  (12)  from  Section  3.2.3  and  the  angle  uncertainty 
€ret  = 0.0235.  The  uncertainty  for  each  retardance  measurement  is  listed  in  Table  9. 


Table  9.  Retardance  error  due  to  retarder  misalignment. 


SRI 

SR2 

SR3 

SR4 

SR5 

89.2 

89.41 

89.4 

87.89 

89.35 

*^RET 

4 X 10-^ 

4 X 10'^ 

4 X 10-^ 

4 X lO'^ 

4 X 10-^ 

The  resulting  error  is  a bias,  and  causes  the  measured  retardance  to  be  slightly  less  than  the  true 
value.  (As  discussed  in  Section  3.2.2,  treating  compensator  error  components  individually  can 
understate  the  total  uncertainty  unless  the  component  uncertainties  are  summed.  Clearly,  the 
uncertainty  is  negligible  for  our  case,  and  our  separation  of  the  errors  poses  no  problem.) 

4.3.2  Retardance  uncertainty  due  to  compensator  retardance  uncertainty.  This  uncertainty  depends  on 
the  specimen  retardance  and  the  compensator  retardance  uncertainty.  We  estimate  a maximum 
compensator  uncertainty  of  3.5°  (Section  2.3.4),  so  application  of  the  error  equation  in  Section 
3.2.1  yields  a maximum  uncertainty  Ad^R.  The  maximum  uncertainty  and  standard  uncertainty  u^-r 
= A6cr/3  is  tabulated  in  Table  10. 


Table  10.  Retardance  uncertainty  due  to  compensator  retardance  uncertainty. 


SRI 

SR2 

SR3 

SR4 

SR5 

AScr 

1.49  X 10'^ 

1.1  X 10'^ 

1.12  X 10-^ 

3.92  X 10'^ 

1.21  X 10-^ 

^CR 

4.97  X 10'^ 

3.67  X 10  " 

3.73  X 10  " 

1.31  X 10-3 

4.03  X 10-" 

The  uncertainty  is  a bias  that  causes  the  measured  retardance  to  be  slightly  less  than  the  true  value. 

4.3.3.  Retardance  uncertainty  due  to  compensator  misalignment.  The  compensator  misalignment  €4 

directly  alters  the  retardance  measurement  with  standard  uncertainty  Uc^,  found  using  eq  (10)  and 
6^,.  Interestingly,  numerical  evaluation  of  the  error  yields  Ucai  = 264.  Table  1 1 lists  the  resulting 
uncertainties. 

4.3.4  Retardance  uncertainty  due  to  Aa  determination  uncertainty.  The  retardance  is  measured  by  direct 
measurement  of  the  magnitude  of  polarizer  rotation  Aa  (Sections  2.2.5  and  2.2.7).  Since  6^, 


11 


= 2Aa,  the  standard  uncertainty  is  UpR,  = 2e^^.  Table  1 1 lists  the  standard  uncertainty  for  each 
measurement. 

4.3.5  Retardance  uncertainty  due  to  compensator  misalignment  (second  measurement^).  Uncertainty  is 
calculated  as  in  Section  4.3.3;  again  numerical  evaluation  shows  that  the  compensator 
misalignment  directly  alters  the  retardance  measurement.  For  the  second  measurement,  additional 
uncertainty  accrues  (Section  4.2.6),  and  the  resulting  standard  retardance  uncertainties  Uca2  = 2eg 
are  hsted  in  Table  11. 

4.3.6  Retardance  uncertainties  due  to  Aa  determination  uncertainty  (second  measurement!.  The 
retardance  is  again  measured  by  direct  measurement  of  the  magnitude  of  polarizer  rotation  Aa 
(Section  4.2.7)  and  the  standard  uncertainty  is  UpR2  = 2e^„.  (Table  11). 


Table  11.  Retardance  uncertainty  arising  from  Aa  error. 


SRI 

SR2 

— ... — ^ 

SR3 

SR4 

SR5 

UcAI 

0.0370 

0.0370 

0.0370 

0.0365 

0.0252 

UpRl 

0.0512 

0.0408 

0.0588 

0.0333 

0.0411 

UcA2 

0.0469 

0.0469 

0.0469 

0.0465 

0.0382 

UpR2 

0.0532 

0.0395 

0.0588 

0.0333 

0.0411 

5.  Analysis  of  Rhomb  Retardance  Uncertainty 

In  addition  to  uncertainty  introduced  by  our  null  measurement  instrument,  there  are  retardance 
uncertainties  that  arise  from  the  retarder  itself.  Our  stable  rhomb,  while  more  nearly  achromatic  and  angle 
insensitive  than  typical  retarders,  exhibits  small  retardance  changes  with  wavelength  and  incident  angle. 
An  estimate  of  these  changes  is  important  because  these  devices  were  blindly  measured  by  each  of  three 
independent  methods  to  ensure  the  reproducibility  and  accuracy  of  our  measurement  techniques. 
Uncertainties  due  to  rhomb  retardance  changes  are  estimated  below. 


5.1  Retardance  Uncertainty  Due  to  Input  Angle 


The  angle  of  total  internal  reflection  (TIR)  inside  the  rhomb  is  0,  = ap  + sin  '(sin((3)/n),  where  ap  is  the 
acute  angle  of  the  rhomb,  P is  the  angle  deviation  from  normal  incidence  at  the  input  face,  and  n is  the 
refractive  index  of  the  rhomb.  Our  device  is  made  with  two  rhombs;  for  the  two  reflections  in  the  first 
rhomb  the  TIR  angle  is  Oj,  while  the  second  two  reflections  in  the  second  rhomb  have  TER  angle  02  = cCr 
- sin  ’(sin(P)/n).  The  retardance  R upon  TIR  is  given  by  Fresnel’s  equations  and  can  be  written  as 


/ 

/?(0)  = 2tan'^ 

\ 


\ 

y/ sin^(6)  - Ifn  ^cos(0) 

sin^(0)  ^ 


(16) 


The  total  rhomb  retardance  is  then  2R(0i)  + 2R(02),  and  the  retardance  error  due  to  input  angle  error  p is 
A6  = 2R(0])  + 2R(02)  - 4R(aR). 
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At  1.32  jim,  n ==1.805  for  the  glass  used  in  out  rhombs  [4].  We  estimate  the  maximum  P by  assuming  that 
the  rhomb  can  be  aligned  using  back-reflection  of  an  incident  beam  with  a resolution  of  one  beam 
diameter.  For  a nominal  path  of  50  cm,  and  a maximum  beam  diameter  of  3 mm,  the  maximum  input 
angle  uncertainty  P^j^^  = tan' *(3/500)  = 0.006  rad.  The  acute  angles  at  which  each  rhomb  was  cut  and 
corresponding  uncertainty  estimates  Up  are  listed  in  Table  12.  The  standard  uncertainty  is  one-third  the 
maximum  uncertainty  resulting  from 


Table  12.  Acute  angles  and  input  angle  uncertainty  for  measured  rhombs. 


Device 

angle  Kr  (°) 

A6(°) 

Up(°) 

SRI 

76.375 

SR2,  SR3 

76.391 

4.17  X 10-^ 

1.4  X 10'" 

SR4,  SR5 

76.366 

SR6 

76.295 

4.21  X 10" 

1.4x  10-" 

Positive  and  negative  values  of  input  angle  P both  decrease  the  retardance  of  the  rhomb. 

5.2  Retardance  Uncertainty  Due  to  Wavelength  Changes 

Using  the  Sellmeier  equation  for  refractive  index  [3]  and  eq  (16),  we  can  estimate  the  rhomb  retardance  at 
various  wavelengths.  For  a maximum  wavelength  variation  of  ±20  nm  about  a center  wavelength  of  13 10 
nm,  the  maximum  retardance  change  is  ± 0.0105°  and  the  standard  uncertainty  is  U;^  = 0.0035°. 

5.3  Retardance  uncertainty  due  to  error  in  rhomb  alignment 

The  stable  retarder  is  made  by  attaching  two  rhombs  with  retardance  During  the  alignment  it  is 
possible  to  include  a small  rotation,  or  twist,  €(  between  the  devices.  Twist  causes  an  error  in  the  retarder 
alignment  €r£j;i-  = - V2e^  since  the  nulhng  angle  between  crossed  polarizers  is  found  with  less  resolution; 
this  value  is  generally  different  than  the  error  €^31-  = 0.0235°  due  to  polarizer  extinction  and  rotation  errors 
(Section  4.2.3).  The  resulting  retardance  error  Urjj  is  calculated  from  the  equations  in  Section  3.2.3  as 
before  and  tabulated  in  Table  13. 


Table  13.  Uncertainty  due  to  twist  in  rhomb. 


SRI 

SR3 

SR4 

SR2 

SR5 

e, 

0.36° 

0.6° 

0.72° 

0.24° 

0 

d 

^RETn- 

0.18° 

0.3° 

0.36° 

0.12° 

0.1° 

Uret 

2.2  X 10'^° 

6.3  X lO  "' 

8.9  X 10'^“ 

1.0  X 10-^° 

6.9  X 10-"° 

These  values  are  larger  than  the  Uj^  values  calculated  from  polarizer  errors  in  Section  4.2.3;  we  will  use 
these  larger  values  in  the  combined  uncertainty  calculations  below. 

6.  Net  retardance  uncertainty  from  measurements  using  null  method 

We  now  combine  the  various  retardance  uncertainties  previously  discussed. 


13 


6.1  Retardance  Measurement  Uncertainty  Arising  from  Effects  that  Bias  Retardance 

Retarder  misalignment  (Uret,  Section  4.3.1)  and  compensator  retardance  uncertainty  (u^r.  Section  4.3.2) 
cause  measurement  biases  that  result  in  a retardance  estimate  that  is  slightly  less  than  the  true  value. 
Similarly,  input  angle  uncertainty  (Up,  Section  5.1)  also  decreases  retardance.  These  are  combined  using 
the  root-sum-of-squares  (RSS)  method  to  give  the  combined  uncertainty  Uc/s  arising  from  biasing.  Table 
14  lists  these  values. 


Table  14.  Uncertainty  arising  from  systematic  effects. 


SRI 

e, 

SR2 

SR3 

SR4 

SR5 

Uret 

2.2  X 10-^ 

1.0  X 10'^ 

6.3  X 10'^ 

8.9  X 10-^ 

6.9  X 10-" 

^CR 

4.97  X 10-' 

3.67  X 10  " 

3.73  X 10-" 

1.31  X 10'^ 

4.03  X 10-" 

Up 

1.4  X lO"* 

1.4  X 10-" 

1.4  X 10" 

1.4  X 10-" 

1.4  X 10-" 

*^c/s 

2.6  X 10'^ 

1.1  X 10-^ 

6.3  X 10-^ 

9.0  X 10-^ 

8.1  X 10-" 

6.2  Retardance  Measurement  Uncertainty  Arising  from  Random  Effects 

Random  uncertainties  arise  from  compensator  misalignment  (Ucm  and  u^^.  Sections  4.3.3  and  4.3.5),  Aa 
determination  (upRj  and  UpR2,  Sections  4.3.4  and  4.3.6),  and  wavelength  variation  (u^.  Section  5.2).  These 
are  tabulated  in  Table  15,  and  the  combined  uncertainties  due  to  random  effects  for  each  measurement, 

Uc/R,  and  u^/r^,  are  calculated  using  RSS.  The  measure  retardance  8^,  is  the  mean  of  these  two  measured 
values,  and  the  combined  uncertainty  for  these  two  measurements  gives  the  combined  standard  uncertainty 
from  random  effects,  Uc/r  = + Uc/r2^)'^-  The  results  are  listed  in  Table  15. 

Table  15 

. Uncertainty  arising  from  random  effects. 

SRI 

SR2 

SR3 

SR4 

SR4 

UCAI 

0.0370 

0.0370 

0.0370 

0.0365 

0.0252 

UpRl 

0.0512 

0.0408 

0.0588 

0.0333 

0.0411 

0.0035 

0.0035 

0.0035 

0.0035 

0.0035 

^C/Rl 

0.0633 

0.0552 

.0695 

0.0495 

0.0483 

UcA2 

0.0469 

0.0469 

0.0469 

0.0465 

0.0383 

UpR2 

0.0532 

0.0395 

0.0588 

0.0333 

0.0411 

Uc/R2 

0.0709 

0.0614 

0.0752 

0.0572 

0.0561 

0.0475 

0,0413 

0.0512 

0.0379 

0.0371 
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6.3  Reported  Uncertainty  for  Null  Measurements 

Comparison  of  Tables  14  and  15  shows  that  systematic  uncertainty  u^/s  is  negligible  compared  to  the 
standard  uncertainty  u^/r  arising  from  random  effects  for  devices  SRI,  SR2,  and  SR5,  so  correction  of  the 
retardance  measurement  using  u^/s  will  not  be  attempted.  For  SR3  and  SR4  the  value  u^/j  be  added  to 
the  measured  retardance  value  as  a correction.  The  combined  standard  uncertainty  Uc  for  these 
measurements  is  u^/r.  Because  u^/r  is  much  larger  that  the  correction  Uc/s  added  to  the  retardance  values  of 
SR3  and  SR4,  we  consider  the  uncertainty  from  the  correction  to  be  negligible.  In  Table  16  the  results  of 
our  prototype  SRM  retardance  measurements  using  the  null  method  are  listed.  In  this  Table,  the  number 
following  the  ± symbol  is  the  expanded  uncertainty  U = kuc  where  a coverage  factor  k = 2 is  used.  If  the 
uncertainty  is  approximately  normally  distributed  with  approximate  standard  deviation  u^,  the  value  of 
retardance  lies  within  the  interval  defined  by  U with  a confidence  level  of  approximately  95%  [6]. 


Table  16.  Results  of  retardance  measurements  using  null  method. 


SRI 

SR2 

SR3 

SR4 

SR5 

Retardance  ( ° ) 89.2  ± 0. 1 

89.41  ±0.08 

89.41  ±0.10 

87.90  ± 0.08 

89.35  ±0.07 

7.  Conclusions 

Measurements  of  a stable  retarder  using  a manual  null  polarimeter  have  an  expanded  uncertainty  U < 0.1  °. 
The  system  is  fairly  easy  to  align  and  operate,  though  two  components  are  wavelength  sensitive  and  must 
be  changed  if  measurements  at  other  wavelengths  are  required.  Angle  readout  limitations  of  the  rotation 
stages  used  in  our  system  contribute  the  largest  error  sources,  and  this  could  be  mitigated  by  improved 
stage  resolution  through  automation  (encoded  stepper-motor  driven  stages,  for  example). 

The  retardance  measurements  used  in  this  uncertainty  analysis  were  performed  by  Ian  Clarke  (now  at 
Galileo  Electro-optics,  Sturbridge,  MA)  during  his  visit  as  a Guest  Scientist.  We  gratefully  acknowledge 
his  work  on  his  project. 

8.  References 

[1]  Tompkins,  H.G.  A user’s  guide  to  ellipsometry.  Academic,  Boston  : 26-30;  1993. 

[2]  Azzam,  R.M.;  Bashara,  N.M.  Elhpsometry  and  polarized  light.  North-Holland,  Amsterdam:  166-169; 
1977. 

[3]  Rochford,  K.B.;  Williams,  P.A.;  Rose,  A.H.;  Clarke,  I.G.;  Hale,  P.D.;  Day,  G.W.  Standard 
polarization  components:  progress  toward  an  optical  retardance  standard,  in  Polarization  analysis  and 
measurement  11 , Proc.  SPIE  2265:  2-8;  1994. 

[4]  Rochford,  K.B.;  Rose,  A.H.;  Williams,  P.W.;  Wang,  C.M.;  Clarke,  I.G.;  Hale,  P.D.;  Day,  G.W. 

Design  and  performance  of  a stable  linear  retarder.  Submitted  to  Applied  Optics. 

[5]  For  example,  see  Kliger,  D.S.,  Lewis,  J.W.;  Randall,  C.E.  Polarized  light  in  optics  and  spectroscopy. 
Academic,  Boston;  1990. 

[6]  Taylor,  B.N.;  Kuyatt,  C.E.  Guidelines  for  evaluating  and  expressing  the  uncertainty  of  NIST 
measurement  results.  Natl.  Inst.  Stand.  Technol.  Tech.  Note  1297;  1994. 


15 


'' V'.  .vd-i,U  ^ 

■ >. ..  , , ’t' 

Va  «» fiU.i'^'s/  ^r’J  Ql  > ■■■  ,u  T<sytal'  «»' 

'..  ili/«i> ‘..iri’dt  ^WiiTiii  ^ '-rfi 

^yijMun  sdJ  .:i!k|,f;''r  4^1;  £j  I?  HK>. 

'<;;i  U ii  tr?f  baJNCnsq/.Jik  ii ''(titt^  ± ©rf^vi'^oUol 

, . ^0.  l'-<v9l  ,*’.M  (w«3t8£0 ii  rttii(|.\J^(lihUtt^  wiynirl^jw  *4iX‘‘J®rtiSto^£ii!!>i 

k t)..!  ><  ilO* 

'.^'r:,  'la-  ^ tjtt?  ■■i^.,g‘  -tmfei'i'  '^^^lar^, 

■ (.-’'s-'j'?  V-  • ^ ti...ii^.'?s-  '^*  ,1V  VyikmtW^^'. 

— • -.  . '..«.w  ...  *»<'g*'*  ■■■. » IIK  > ’.'^  ■**—....  .1.1..  . — ■Ite.i.i^.r  ' ~~  ■ »1W1»||..  <.ii  I j-  - ^1  . |»««.|  !■•  nlur  Bp  in  i i/MtttH 


I 


) ll>  . 


.**■'■(«.  is-'' 


*.fto}2ublM>^  X 


! 0 i VI  !A<(>n»i*.3[i\  j '.N  S'/^id  T9’'  ncmfcy.  !ii'‘{  (?!W(.ffltyj  • ?t)iw  ftkteJz  6 1o  iJi»fnmwaK»M 

•■r.  ;<►  ■ i’tt,!  . ■ ,■  ;.UV».  . ’-.'-j.iT?  {»,  v/J04)> 


I il  > 


■ *.%■  ,ifT^jte.a^^4  '•i*' 


'»  '.  'Jl- 


I ■■.  r-1  f>'u’S|||f!n»  sH-r4ioc4>.flE, »-  .-.iTir'-vv4ii«i»-  ^ 

{ -itfirCi^  »t>:  f>.  l(  it\,i<8n«:;‘-  - j-’  i '■• , !”^{iir/«jr4i<>ijW0»i»ru»46J« 

,,:«  . » ; " ; . r1'»J;,t  i.twatTtriirci  5r.tAtd(»rii 

• I , Q ivtl  ^ v 'fU  •«#irm-'.5«  rl#*loir‘  ItaA#  '.f»>r-<'«W*  ‘'.  ‘ *’'T»  «on«tm5si  wTf 

. • ; •/■  ■■  I .>•  ‘ ■'l'■^*;J  ■:.'’  t 'f^  Ik-uvrirj'*.  ruisO  > ,y,  ySfDii'idiWft 

■•  ''  ' ■ ■ -JS!.*.  ■ ’« 

— - ...  ...  ;'  BfnJnSffF  M 

,r.  I r?i"  ■ - -'Vsfri^'i  . 

•v''i:l-  U^  ;\i.  fif  .j«W.r-tis:iA  of  iU  W IM 

,a’’i  f-Hut  \ .■'ii’'''itlli 'Hhv)/  ; MJS  %\Mk  m 

tyui^nu<:  . ^ „ •}  ^ ;..ifel‘  ,j<.  '"■•i  , . *> 

' ■ ■^'’ i . r-i  4' 

.r  L’  .veQ  :.da'^%fii%'i.D4  .gnsW  (f‘l 


^t.-.  -■  "Wj; 


■V  •:*  > '-Wv^i’^A  - i‘5tr*w‘'-  - 


